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Abstract 

We study the properties of cold asymmetric nuclear matter at high 
density, applying the quark meson coupling model with excluded vol- 
ume corrections in the framework of the Landau theory of relativistic 
Fermi liquids. We discuss the role of the finite spatial extension of 
composite baryons on dynamical and statistical properties such as the 
Landau parameters, the compressibility, and the symmetry energy. 
We have also calculated the low lying collective eigenfrequencies aris- 
ing from the collisionless quasiparticle transport equation, considering 
both unstable and stable modes. An overall analysis of the excluded 
volume correlations on the collective properties is performed. 

PACS : 12.39.Ba, 21.30.Fe, 21.65.-Ff, 71.10.Ay 
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1 Introduction 



The study of the nuclear medium composed of different fractions of protons 
and neutrons has been developed for long time and it has concentrated a 
renewed interest in the last years. The equation of state of isospin asymmet- 
ric nuclear matter is a subject of particularly intense research [H [21 [Sj HI [5] . 
The possible applications range from the structure of radiative nuclei, the 
dynamics of rare isotopes and the cooling process of neutron stars. A key 
role in many of these calculations is played by the density dependence of the 
symmetry energy [3j, which may be extracted from recent isospin diffusion 
data in heavy ion collisions experiences [2j. 

The standard theoretical calculations have used as the relevant degrees of 
freedom protons and neutrons, moving non-relativistically through poten- 
tials representing the averaged instantaneous interaction among nucleons. 
Furthermore, since the development of the Quantum Hadro-Dynamics f6], 
it has became a common practice to use relativistic hadronic fields in these 
kind of calculations. In such the nuclear interaction is mediated by 

mesons of scalar or vector isospin character. The role of the scalar isovector 
meson (oq (980)) has been emphasized in recent investigations [7J. 

A further step has been given by models incorporating the quark struc- 
ture of hadrons. They may be reduced to mean field hadronic pictures whose 
parameters, such as masses and vertices, hide the quark dynamics. These ef- 
fective models provide both a connection between the hadronic phenomenol- 
ogy and the fundamental theory of the strong interactions, besides a sound 
description of hadronic matter and atomic nuclei. 

Among the dual quark-hadron theoretical frameworks, those based on the 
bag models include explicitly the quark confinement volume. Specifically, 
the Quark Meson Coupling (QMC) [8J states the dynamical evolution of 
this confining region, which depends on the global properties of the nuclear 
medium as well as on the configuration of the hadronic fields. However, this 
feature is usually missed [8j in passing to a pure hadronic context, as baryons 
and mesons are regarded as point-like particles in this limit. 
The relevance of the finite extension of nucleons in the evaluation of some 
nuclear statistical properties has been stressed long time ago [9] . In order to 
fulfill this assertion, several corrections have been introduced in the hadronic 
interactions, mainly invoking a Van der Waals-like normalization pUl [TT| 
[T2l [13] . The motivation for such a procedure is to parameterize in compact 
form the strong baryon-baryon repulsion at very short distances. In the high 
density realm of nuclear matter the mean free path of nucleons is comparable 
to the quark confining size, hence excluded volume effects become relevant. 
The correction due to the finite extension of nucleons have been introduced 
into the hadronic sector of the bag models [HI [13] , it was found that it is 
necessary to describe properly high density hadronic matter without violating 
the model assumptions jT3] . 

Another theoretical scheme suitable to describe the nuclear collective phe- 
nomena is the Landau theory of Fermi liquids. Although it was originally 
stated in a non-relativistic fashion, see for instance [H], it was subsequently 
extended to deal with the relativistic fields formalism flbi [T6] . The Landau 



parameters are useful to evaluate thermodynamical properties, the stabil- 
ity conditions in phase transitions, the nuclear matter collective excitations 
which couple, for instance, to the weak interaction governing the neutrino 
emission of URCA processes in neutron stars. 

In this work the behaviour of isospin asymmetric nuclear matter is studied 
in the framework of the QMC model with excluded volume corrections in 
the hadronic sector. Special attention is paid to the Landau parameters and 
the collective modes. In the next section we give a resume of the QMC 
model and we describe the correlations generated by the Van der Waals-like 
normalization. In section 3 we define the Landau parameters and use them 
in section 4 to derive the isoscalar and isovector collective modes. Results 
and discussions are presented in section 5, and finally the conclusions are 
drawn in section 6. 

2 The Quark Meson Coupling Model 

QMC is an effective model of the quark structure of hadrons, it is inspired in 
the MIT bag model, so that in its starting point the confinement mechanism 
has being accomplished and chiral symmetry has been broken. In order to 
describe the dynamics of the emerging hadrons, usually the model is pro- 
jected into a picture of point-like baryons interacting through virtual mesons 
[8j. the same ones which couple quarks inside the confinement region. Since 
wc are primarily concerned with isospin asymmetric nuclear matter, we only 
consider u and d flavor of quarks coupled by the isovector p and 6 (ao (980)) 
mesons, besides the commonly used iso-scalar a and u ones. 
Within the QMC model baryons are represented as non-overlapping spheri- 
cal bags containing three valence quarks; the bag radius changes dynamically 
with the fields configuration. 

The mean field approximation (MFA) is the natural scheme of solution, which 
replaces the meson fields by its classical expectation values. These mean val- 
ues determine the nucleon fields, which in turn become the source of the 
meson ones. 

In the MFA the Dirac equation for a quark of flavor q, [q = u,d), of current 
mass rriq and Iq third isospin component, is given by [8] 

- gl^'u - g^Iq - mq*)^^ = 0, (2.1) 
here the notation rriq* = niq — g^a — Iqg'^S is used. 

For a spherically symmetric bag of radius Rb corresponding to a baryon of 
class b, the normalized quark wave function \l/^(r, t) is given by 

where x'' is the quark spinor and 

e,b=^ + glco + gllq p, (2.3) 



= Rb' [2ng,{n,, - 1) + R,mg*] (2.4) 
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^qb — Rb^riq* 
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with VLqb = [x'^^ + (i?bmq*)^]^/^. The eigenvalue Xqb is the lowest solution of 



joi^qb) ^ PqbjliXqb), (2.6) 

which arises from the boundary condition at the bag surface. 

In this model the ground state bag energy is identified with the baryonic 

mass Mb, 

M,^5i!!5|iZi?! + W_ (2.7) 

Kb o 

where is the number of quarks of flavor q inside the bag. The bag constant 
Bq is numerically adjusted to get deflnite values for the proton bag radius, 
and the zero-point motion parameters zob are fixed to reproduce the baryon 

spectrum at zero density. 

The dispersion relation for the b-baryon is 

k'o = ^M62 + (k-E,)2 ± g^uo ± Qphpl, (2.8) 

for particle (+) and antiparticle (— ) solutions. Within the MFA at zero 
temperature only the particle solutions contribute. In this expression we have 
assumed that the strength of the couplings does not depend on the quark 
flavor. We have also introduced the baryonic isospin projection 7^ = X]q '^q^?; 
and the vector nucleon sclfcnergy = Qi^lj + Qpph- 

For homogeneous infinite static matter the spatial dependence of all the 
meson fields can be neglected, so that its equations of motion reduce in the 
MFA to 



Y V d(T Jr,s 



5 = 



^E^-^.'' (2-10) 



^2Y.9ph3i. (2.12) 



Here = (n , j ) stands for the mean value of the nucleon current density 
of isospin b. In the reference frame where the averaged momentum of matter 
is zero the spatial components of the currents become null j = 0, hence only 
the time-like projections of the vector mesons are non-zero. 



It must be noted that Eqs. (12. lip and (I2.12p refer only to the third isospin 
component, since the remaining ones become zero in the MFA. 

From the relations (12.71) and (12.80 it can be seen, as it was earlier men- 
tioned, that the nucleon mass and energy spectrum depends on the assumed 
meson classical values. They are in turn, determined by the nucleon densities 
as given by the equations (I2.9p - (l2.12p . 

The densities and j^^ are defined with respect to the ground state of the 
hadronic matter, which at zero temperature is composed of baryons filling 
the Fermi sea up to the state with momentum kpb 

= ^ e{kF, - |k|), (2.13) 



In Eqs. (I2.13P and (12.141 ) the factor {} is included for future use and it takes 
the value = 1 for point-like baryons. 

The total energy density S and pressure Pq of hadronic matter for point- 
like baryons is evaluated as 

rkpt 



+ — y / dkkWM,^ + k^, (2.15) 
Po = (2.16) 



where /Iq = kQ^kph), see Eq. (12. 8p . is the chemical potential for point-like 
baryons. 

In the QMC the radius Rb is a variable dynamically adjusted to reach the 
equilibrium of the bag in the dense hadronic medium. We use the equilibrium 
condition proposed in ref. [13], which can be obtained by minimizing the 
energy density £ with respect to Rb 

_ J_ fdMb] =J_^ dkk' 12 m 

AnRl \ dRb ) 37r2 ^ ^ J, ^M,' + k^' ^ ' ^ 

where ^ = 1. This result refiects the balance of the internal pressure of the 
bag with the baryonic contribution to the total external pressure, represented 
by the left and right sides of Eq. (I2.17P , respectively. 

The factor ^ will be redefined below where excluded volume effects will be 
considered. 

The QMC model heavily relies on the assumption of non-overlapping 
bags, using this criterion an upper density limit around three times the sat- 
uration density of symmetric nuclear matter has been found [T7|. In order 
to properly take into account this severe restriction, a simplified model was 
introduced in [13], which describe baryons as extended objects. Since finite 



size baryons are assumed non-overlapping, their motion must be restricted 
to the available space V defined as [TO] 



V' = V -"^N\, (2.18) 

b 

with A^* the total number of baryons of class b inside the volume V, and 
the effective volume per baryon of this class. The last mentioned quantity is 
proportional to the actual baryon volume, i.e. for spherical volumes of radius 

Rb 

4:71 

Vh = a—Rt', (2.19) 

where a is a real number ranging from 4, in the low density limit, to 3\/2/7r, 
which corresponds to the maximum density allowed for non overlapping 
spheres, in a face centered cubic arrange. Since we wish to study the high 
density regime of homogeneous isotropic matter, we shall adopt a = 3\/2/7i 
in all our calculations. 

Consequently it was assumed that the nucleon fields can be normalized 
replacing V for V, which in turn implies a normalization of the nucleon 
densities f l2.13p . (12.141) and of the nucleon contribution to the energy (I2.15p . 
The final result may be obtained by taking = 1 — 'Y^^n^Vb within these 
equations and ^ = in!]./{aii!') in Eq. ( 12.17^ [T3j . 

Since d introduces an explicit dependence upon the baryonic densities, 
the chemical potentials get an extra term, i.e. 

/ = =/iUA;z^ (2.20) 



dn, 
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Vb ^ f''^^ dkk^ 



A^^ = , ■ (2.21) 

Correspondingly, the total pressure acquires an additional term AP as 
compared to the pressure of point-like baryons Pq in Eq. (12.161) 

Ph = Po + AP = Po + ^ n^A/. (2.22) 



3 Landau parameters and equation of state 

The Fermi liquid theory of Landau assumes that the low-lying excitations 
of a physical system admits a representation in terms of quasi-particles and, 
circumstantially, collective modes. If the quasi-particle states can be iden- 
tified by a composed label B = (6, /?), where we have singled out the first 
place for the isospin projection and (3 collects discrete spin and momentum 
indices, then the occupation number of such a level is denoted by fs- The 
conserved baryonic number can be expressed in terms of a summation over 
such distribution functions: n = '^^fs- The same statement is valid for 
every extensive conserved quantity such as the energy, which can be written 



as J2b fs^B plus current-current interactions, where is the nucleon sin- 
gle particle spectrum. For infinite nuclear matter the summation over the 
discrete momentum indices must be replaced by an integration over the con- 
tinuous spectrum with measure d^k/{2TT^). 

The formulae fl2.13l) and fl2.14p for the nucleon densities must be rewritten in 
order to fit this form. Firstly the replacement Pi = ki 'd^/'^ is made, secondly 
we identify Y.spin ! d^^pph - |p|)/ (27r)^ J2pfbf3, with ppb = kp^^^^^- 
Hence, Eqs. f l2.13p . (12.141) . and (12.151) can be rewritten as 



^ yM^ + (^-i/3p^-S,)2' 
V VM2 + (^9-V3p^-S,)2' 

B ^ ' 

+ 1 (^2 ^2 ^ ^2 ^2 _ ^2 ^^^^ _ ^2 ^^^M) ^ (3.3) 

where the momentum independence of the nucleon masses, characteristic of 
the Hartree approximation, has been emphasized. 

In the last formula E must be considered as a function of the distribution 
functions /s, the bag radii Rb and the meson fields = a, tu, 5, p, hence first 
variation of the energy density yields 

i?, 0) = E ^^^^ + E i^'^ + E 1^ ^z^- 

^ dRb ^ B 

The first term is null because of the equilibrium condition for the bag in the 
nuclear medium, Eq. (12.171) . the second one is hkewise zero due to the meson 
field equations (I2.9p - (l2.12p . Therefore, the quasiparticle energy spectrum 
according to the Fermi hquid prescription is 

d£ I 

£b = jr^ = jM^ + {i3-y^pp-i:by + g^Uo + gpPoIb 

OjB ^ 

Vb^ ^-'I'VP> ■ {^-"'vp> - S.) 

If Eq. (13.41) is evaluated at the Fermi surface, i.e. |p/3| = pp = PFb, in the 
limit of isotropic matter and restoring the continuum spectrum, it yields the 
chemical potential given by Eqs. (12.201) . (12.211) . in support of the thermody- 
namical consistency of our approach. 

In order to evaluate the Landau's parameters of the nuclear interaction, 
a second variation must be performed, i.e. Fb' b = d'^E /dfB' dfB = Fbb'- 
Since all of the bag radii and the meson fields depend ultimately on the dis- 
tribution functions /b, Eqs. (12.90 - 02.120 and (12.170 must be differentiated in 
order to obtain a closed expression. 



The Landau's parameters are defined as the Fourier coefficients of an expan- 
sion in terms of the Legendre polynomials 



F'sB' = (/ + 1/2) j dxPi{x) 



where the integrand must be evaluated on the Fermi surface at the end of 
the calculations. 

It is found that the second order derivative in the expression above, depends 
linearly on the variable x = p/3/ ■ Pp/iPis' P/s), hence all the parameters of 
order greater than one are null, whereas for / = 0, 1 we have 

c=n,p c \ a sc / 



9^ 

q 



+^'^^^^ + E (<Qi^^ - y E <QqcMcHP ) y,,, (3.5) 

q=u,d \ ° c=n,p , 



pi ^ Kh + lahCp + 2C^Cp{l + Iah)na/Ea 

EaE, 1 + (a + Cp) Ec=p,n t + 4C.C,ff 



(3.6) 



Here the definitions = {g^/m^y for (f) = uj, p and Ea = a/ + k\^, have 
been used, together with 

^ / ^(3) _ f^Fa 



7r2(M2 + P)'/2' » 7^2^, 

qa ^ ) 

^; 5 a 



(jra 



Furthermore, the variables Xah and Yqi, are the derivatives = dRa/Ofs 
and = drrtq* /dfB evaluated in the isotropic matter limit and taking all 
the momenta at the Fermi surface. Actually, they are the solutions of an al- 
gebraic coupled system of equations obtained by differentiating Eqs. fl2.17p . 
f l2.9l) . and (12. lip and making a linear combination of the two last results. 
For further development, we define symmetric adimensional Landau param- 
eters 

•^afe = ^ba = Vr^TfoF^^ (3.7) 

where F;, = (kphEb/n^) (6 = p,n) is the relativistic quasiparticle density of 
states at the Fermi surface. 

Eqs. (I3.5p and (13.60 coincide with the previous results [151 [16] ^oi symmetric 
nuclear matter if the limit —>■ is taken carefully. 

The Landau parameters can be related with the nuclear compressibility 
K and symmetry energy Es- For a given nucleon density n = + and 
asymmetry coefficient t = (n" — n^) / n, it can be shown that 



K, = 9n 



9n 



ty I (i+tr 



and 



2n\dt 



n 



8 vr 



r 



(3.8) 



(3.9) 



where the effective Landau parameters and are, respectively 



(1 - tf r„ + 2(1- 1^) yr;r;^°„ + (i + tf v, j^l 



p TtO _ 9 /p p -ttO _i_ p -ttO 

* ~ r + r 



(3.10) 



The symmetry energy Es{n) is obtained from Eq. fl3.9p in the hmit of t = 0. 
These expressions reduce to the usual ones in the symmetric case without 
volume corrections {d = 1), see for instance [TB] . 



4 Collective excitations 

Collective modes are associated to local density fluctuations that propagate 
in the hadronic mean field. These fluctuations are the effect of small pertur- 
bations of the occupation distribution fs around the nucleon Fermi level of 
the quasiparticles. At zero temperature the unperturbed nucleon distribu- 
tions are given by [H] 

fB = ^e{fI,-6B) (4.1) 

where V is the volume of the system, and fib is the chemical potential. Local 
density fluctuations add a small variation of the occupation numbers around 
the fs equilibrium value, i.e. \1A^ il5j 

/B(r,t) = fB + SfBir,t) = fB + u^e'^^-^--'^ (4.2) 

and to first order the quasiparticle energies Eb change accordingly 

6eB{r,t) = Y,FBB'5fB'{r,t) (4.3) 

B' 

The propagation of these perturbations at low temperature is governed by 
the collisionless Landau's kinetic equation 

dfs ^ dfB_ deB_ _ df^ deB_ ^ ^ 
dt dr dpfs dpi3 dr 

Introducing Eq. (14. 2 p into Eq. (14. 4p and keeping only linear contributions of 
the fluctuations, we obtain 

dSfB ^ dSfB dsB _ dfB_ dSsB ^ ^ 
dt dr dpp dpf3 dr 



which can be further reduced to 



Ub 



vr^ ^ — ^ 



(4.6) 



a=p,n 



The sum should be done at the Fermi surface of protons and neutrons. For 
isotropic matter the direction of propagation can be arbitrarily choosen along 
the azymuthal axis, so that (pF6-q) = PFb <? cosxb- Expanding the ampli- 
tudes Ub and the Landau parameters Fba in terms of Legendre polynomials 



Ub 



^ uiPi (cosxfe) 



-^FbaPi{cos{xa-Xb)) (4.7) 
and making use of the addition theorem for Legendre polynomials, we have 

5^ FfcaWa = 5^ 5^ m^PL<Pl (cOSXfc) (4.8) 



a=p,n 



a=p,n I 



Replacing in Eq. (14. 6 p and writing all back in terms of the Fermi momenta 
kpbi we arrive to the following system of homogeneus linear equations 



{21 + 



a=p,n 



(4.9) 



with Sb = ^^^"^ {uj/q)/vFb-i where vpb = (kph/Eb) is the relativistic quasiparti- 
cle Fermi velocity. The function Qui is given by 



Qii'{s) = Qm{s) -- 
which take the particular expressions [1 

fioo(s) ■ 



y 



dyPiiy) , . 
1 yy-s) 



Pv{y) 



(4.10) 



1 + - In 

2 V s + 1 



VLii{s) = snio{s) + -6n 



(4.11) 



In the case of unstable modes for which s becomes a purely imaginary 
quantity in the upper half complex plane, i.e. s = i( {( > 0), we have 
^oo(^C) = 1 ~ C arctan(l/^). 

Using Eq. (14. lip and keeping in mind that only terms with / = 0, 1 are non 
vanishing, we can solve Eq. (14. 9 p for the amplitudes Up and m° 



u" 



u^ 



Cp ^ 
o Sp y s 



p \ ^ n 



Y 

-^T^ u^ 

-L r 



Cn 1 ^ 



J^^ U^ 

^ rip 



(4.12) 



where Cb = {1 + | JF^^), (b = p,n). Replacing into the equations for / = 1, 



we have 



A 7,1 -u . /li A 7/1 - n 

V -L 7?, 




-^A„pmJ + A„„m^ = (4.13) 



where 

App = Cp + i)QQQ{Sp){CpJ-'pp + '&- ^pn^np~^ ^p^ ^pp) 

6 Vpp 

\n = i)[\:Fl^ + noo{sp){'^TlpTl^ + '^c^^^^^ 

The remaining coeficients Ann-, Anp are obtained through the replacement 
p ^ n in the previous formulas. The nontrivial eigenmodes u\{m)^ (m = 
1,2, ...) of the linearized transport equation (14.91) are equivalent to the non 
vanishing solutions of Eq. (l4.13l) . Therefore the corresponding eigenvalues 
s^^^ satisfy 

■^pp -^n n -^pn -^np (4. 15) 

For fixed nucleon density n and isospin asymmetry t, the solutions s^^^ of 



m- 



Eq. (I4.15P determine the zero sound dispersion relation {oj/q) for the 
eigenmode. Following [18] we shall use the sign of the relative amplitudes 
Q = Up/u'^ to determine the isoscalar {g > 0) or isovector {g < 0) character 
of each mode, for arbitrary isospin asymmetry. 



5 Results and discussion 

Within the present model the current quark masses have been chosen as 
rriu = rrid = 5 MeV, the bag parameter has been fixed as B^^^ = 210.89 MeV 
neglecting its eventual density dependence. The parameter zob was adjusted 
in order to reproduce at zero density the empirical value of the nucleon mass 
and a nucleon bag radius Rb = 0.6 fm. Numerical values for the meson 
masses have been taken as m„ = 550 MeV, = 783 MeV, = 984 MeV 
and rup = 770 MeV. 

Since mesons interact directly with quarks, the corresponding meson-nucleon 
couplings are related to the quark-meson ones g'^ = {(p = a, u, S, p) in a 
simple way by assuming vector meson dominance, i.e. 

gl = 3g:, gt = 3g:, g's = gl g'p = g;, 

for b = p,n. Their numerical values are obtained by reproducing the sym- 
metric nuclear matter saturation properties, i.e. baryonic density, binding 
energy and symmetry energy 

Ebind = {£/n)^ - Mc^ = -16MeV, 

^^('^o) = ^ f = 31.6MeV, (5.1) 

2no \ dt /i^o 



where M = 938.92MeV/c^ is the average free nucleon rest mass. 
The constraint of the symmetry energy alone is not sufficient to determine 
unambiguously both Qp and gs, instead it establishes a non-linear relation 
between them. In order to restrain their possible variation, we have used 
valuable phenomenological data coming from the analysis of isospin diffusion 
in heavy ion colhsions [3J, as discussed below. 

The values of the couplings are sensitive to either the inclusion or not of the 
excluded volume corrections. Both instances are considered in Table [TJ 
The cases with the excluded volume correction (CC) have been compared 
with calculations without it (NC). It must be emphasized that in the last 
case the values ^ = {} = 1, and A/i^ = A/i" = must be used. 
With this set of parameters the equation of state for asymmetric nuclear 
matter with finite nucleon size has been evaluated. The results are displayed 
in Fig. [T]for the the pressure P and the nuclear compressibility k in terms 
of the baryonic number density for several isospin asymmetries. For sym- 
metric nuclear matter {t = 0) at normal density we obtained k = 340MeV, 
which is higher than the usually assumed values in similar calculations. How- 
ever it should be stressed that recent mass measurements of the pulsar PSR 
J0751+18007 yield M = 2.1 ± 0.2Mq [H], which is hardly compatible with 
the usually adopted value k ~ 250 MeV. Indeed a stiffer equation of state is 
needed to reach this observational constraint. 

In previous investigations the authors found that finite baryonic size are 
responsible for a 12% increment in the maximum mass of a neutron star, ob- 
taining Mjnax = 1.89M0. Therefore the strong baryonic repulsion previous 
to the Quark-Gluon Plasma transition contributes significatively to produce 
a stiffer equation of state, and consequently it would give a phenomenologi- 
cal basis to understand the unexpectedly high value measured for the PSR 
J0751+18007 mass. 

The compressibility becomes negative at low densities leading to thermody- 
namical instabilities, but it grows with increasing t at a fixed density. The 
instabilities disappear for nuclear compositions approaching pure neutron 
matter (t = 1) as a consequence of the repulsive character of the asymmetry 
energy. 

The effect of increasing asymmetry on the binding energy Ehind can be ap- 
preciated in the upper panel of Fig. [2l In particular neutron matter remains 
unbound for all densities. A comparison of the equation of state for the CC 
and NC cases was given in [13] for hadronic matter in /5-equilibrium, the 
conclusions given there can be extended to the case of nuclear matter at 
fixed isospin asymmetry, i. e. density dependence for energy and pressure 
are stiffer in the CC instance. 

In the lower panel of the same figure the asymmetry energy Es is dis- 
played, together with the curves corresponding to the empirical expression 
Es = 31.6 {n/rioy evaluated at the limit values 7 = 0.69 and 7 = 1.05 ob- 
tained in Ref. |2J. It can be seen that our result lies between these curves, 
showing a significative agreement with the 7 = 1.05 case in all the range 
0.5 < n/nQ < 2.5. It must be mentioned that within this model, only a nar- 
row range of values for the pair of couplings gs, gp is able to fit the reference 
value for Esirio) and to produce simultaneously a curve entirely comprised 



between the phenomenological constraints. For very low values of gp it is 
not possible to adjust the symmetry energy at the normal density, increasing 
this coupling yields a stiffer density dependence for Es{n), which quickly goes 
beyond the curve 7 = 1.05 of Fig|21 
Assuming a decomposition 

Es{n) = Es{no) + L{n/no - l)/3 + Ks{n/no - l)Vl8, 

we have found L ~ 95.8 MeV and Kg — 19 MeV. The first quantity agrees 
with the value L = 88 ± 25 MeV found in [Ij, whereas we obtained for the 
combination Kasy = Kg — 6 L = —555 ± 5 MeV in comparison with the 
suggested value Kasy = —500 ± 50 MeV [2J. It must be pointed out that the 
empirical value for L is coherent with the inequality 0.7 < 7 < 1.1, however 
the numerical value extracted from experimental data for Kasy favors a stiffer 
symmetry energy with 1.26 < 7 < 1.3. 

A comparison between NC and CC results yields a enhanced growth for Es{n) 
in the last case, although differences become appreciable for densities higher 
than 2no. For the NC case the values L = 91.4±0.6 MeV and Kasy = -553±2 
MeV have been obtained. 

In the limit of point-like baryons the formulae for the compressibility and 
the symmetry energy reduce to 



which agree with results obtained in relativistic field models with structure- 
less nucleons, with exception of the couplings C^^^ = (g^pQ/m^)'^ (0 = a, 6) 
which includes the factor Q defined below Eq. (13. 6p . Taking into account 
that Q depends on the medium properties through m* and R, density de- 
pendent effective couplings have been obtained, due to the quark structure 
of nucleons proposed in the model. 

The density dependence of the adimensional Landau parameters are plot- 
ted in Figs. [3] (CC) andH] (NC). At sufficient low densities all scalar param- 
eters Tpp, J-'np, and are negative, reflecting the attractive character of 
the effective nucleon interaction. Therefore instabilities in the equation of 
state can appear in this density range. Comparing Figs. [3] and H] we ap- 
preciate that, for densities n/riQ < 1, the general trend of excluded volume 
correlations is to slightly increase the nucleon-nucleon attraction. On the 
other hand in the range n/uQ > 1 the volume corrections enhance the repul- 
sion among nucleons as compared with the respective NC cases. This fact is 
expected since the available volume per nucleon decreases, resulting in a non 
negligible compression of the bags at higher densities [13j . 
The finite volume effects are more evident for the ^F^f,, [b, b' = p, n) compo- 
nents, enhancing their absolute values for all the range of densities, specially 
for n/no > 1. 

Because the adimensional Landau parameters contain density dependent fac- 
tors, J-'p'p and J-'p'j^ decrease for growing isospin asymmetry t as they are a 
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measure of the strength of the in-medium proton interaction. The opposite 
is true for J-'^n- 

The low density hmit of CC results qualitatively agree with others calcula- 
tions [ISl [IS], both for nuclear symmetric matter and for the pure neutron 
case. 

It can be seen that in the sub-nuclear realm of the scalar interaction, the 
in-medium p — n strength overrides the p — p and n — n components. Fur- 
thermore the n — p attractive effect has a wider density range, extending 
beyond the normal value hq. 

Collective quantum fluctuations give rise to proper modes which are so- 
lutions of the eigenvalue equation (14.151) . The corresponding dispersion rela- 
tions (uj/q) are displayed in Figs. [5](CC) and[6](NC), for some typical values 
of the isospin asymmetry t. It must be pointed out that stable modes satisfy 
Sp^n > 1, whereas pure imaginary values correspond to unstable propagation. 
In the last case the quantity \uj/q\ has been plotted. With exception of nearly 
pure neutron matter, unstable modes are always present at very low densities 
and therefore they are practically insensitive to finite volume effects. 

As it was mentioned earlier, collective modes can be classified as isoscalar 
and isovector, according to proton and neutron vibrations being in phase or 
in opposition, respectively. 

Unstable modes are found to be isoscalar in character, reflecting the fact 
that both isospin components simultaneously undergo a liquid-vapor phase 
transition, leading to cluster formation [18]. This associated mechanical in- 
stability is evidenced by the negative sign of the nuclear compressibility in 
this density domain. Since k, comprises density variations at fixed isospin 
asymmetry t, Eq. fl3.8p . it preserves the proton to neutron ratio. 
It is easily verified that for the symmetric case isoscalar modes (stable or 
unstable) also conserve this ratio, namely g equals {dn^ /5n"')^ = [n^/n"') for 
t = 0. This is no longer true for iso-asymmetric matter. 
For growing asymmetry t these unstable modes show a decreasing amplitude, 
because the repulsive n-n interaction progressively dominates, until the in- 
stability vanishes. 

For stable eigensolutions (uj/q) represents the zero sound velocity Vg. In 
the iso-symmetric case there are two stable modes, both in CC and NC in- 
stances. The isovector branch appears at very low densities, whereas the 
isoscalar one starts at n/ng > 1 . For the CC plot they cross each other 
at u/hq ~ 1.51 {n/no ~ 1.9 for NC), but keep their own character because 
isoscalar and mechanical oscillations preserve the same proton to neutron 
ratio for t = 0. Thus, they do not couple to isovector fluctuations, which are 
related to species separation [T8]. 

In general for iso-asymmetric matter there is only one stable branch which 
has a mixed character. In fact, the ratio of proton to neutron amplitudes 
Q changes smoothly from negative (isovector) to positive (isoscalar) as the 
density increases. When the nucleon finite size is considered (CC) the change 
of character takes place at lower densities than in the point-like case (NC), 
as can be appreciated from figures O and El respectively. 
As mentioned before, in asymmetric matter the ratio g for isoscalar modes 
(either stable or not) do not follow the constraint of keeping t constant, as 



mechanical oscillations do \l8l. In fact, the asymmetry of the medium in- 
duces a more complicated scenario where some chemical component is also 
involved, i.e. the relative proton to neutron concentration is modified along 
the quantum mode. The mixed iso-character of the stable modes at finite 
t can be ascribed to this cause. At low densities these collective modes 
are isovector like, indicating a significative species separation. As the den- 
sity grows, pressure (mechanical) effects become dominant and induce the 
change to isoscalar like character. 

For pure neutron matter only one stable mode has been found in this range 
of densities. 

Concerning the velocity of propagation, a rough 30% increment in the CC 
results for Vg are observed in the high density domain as compared to the 
NC ones. As an artifact of the approximation used to solve the Landau's 
kinetic equation, it is found that Vg approaches the velocity of light at den- 
sities n/riQ > 4, in the CC approach. However, for such high densities the 
dissipative effects are non- negligible, resulting in a more involved dynamics 
which is beyond the scope of the present work. 

6 Conclusions 

We have studied the thermodynamical properties and collective modes of 
cold asymmetric nuclear matter. This has been performed within a model of 
structured nucleons, taking into account an appropriate normalization of the 
nucleon fields to prevent overlapping of the quark confining regions at high 
densities [13] . Many body properties of hadronic matter are interpreted in a 
Fermi liquid picture with quasiparticles and collective modes. 
Within this framework the relativistic Landau parameters for iso-asymmetric 
nuclear matter have been evaluated, and an explicit relation with the nuclear 
compressibility k and symmetry energy Eg has been obtained. The overall 
trend of excluded volume correlations at densities n > no is to develop an 
stiffer equation of state, relative to the case where they are absent. 
Different quantities such as density dependence of the symmetry energy, its 
slope, and the asymmetry compressibility Kassy have been evaluated, obtain- 
ing qualitative agreement with empirical estimates. 

It is argued that short range hadronic correlations in the dense medium pre- 
ceding the deconfinement phase transition, which has been parametrized in 
the present work as excluded volume corrections, contribute to understand 
the high neutron star mass measured recently. 

We have also applied the present formalism to the Landau's collisionless 
kinetic equation for small fluctuations around the Fermi level. We have ob- 
tained the eigenvalue equation for instability and zero sound modes within 
this scheme, and have discussed the propagation of these modes in the dense 
hadronic medium. 

The behavior of collective excitations can have important consequences in 
heavy ion reactions, where the formation of fragments is chiefly governed by 
isoscalar fluctuations. 
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9l'' 


u,d 
95 


9^/ 


CC 


5.76314 


2.78280 


5.75150 


4.3500 


NC 


5.99339 


3.00770 


5.12075 


1.5000 



Table 1: The quark- meson couplings used in the case with excluded volume correction 
(CC) and without it (NC). 




Fi gure 1: The pressure of cold asymmetric nuclear matter versus density for the CC case 
(upper panel). The different lines correspond to some typical values of the asymmetry 
parameter t. In the lower panel the corresponding compressibility k is shown for the same 
set of t values. 





Figure 2: In the upper panel the binding energy of cold asymmetric nuclear matter for 
the CC case and for some typical values of the asymmetry parameter t is drawn. In the 
lower panel we plot the asymmetry energy for the CC and NC cases, together with the 
empirical expression Eg — 31.6 {n/noy (7 = 0.69, 1.05) Ref.[2]. 



Figure 3; The adimensional Landau parameters for the CC case. In the left (right) panels 
the scalar (vector) components for the pp, pn, and nn pairs are respectively plotted from 
top to bottom, at several fixed values of t. The line convention valid for all the cases is 
indicated in the central right panel, it must be noted that t = 1 is relevant only for the 
nn amplitude. 




Figure 4: Idem as 



FigEl but 



for the NC case. 




Figure 5: Dispersion relation for collective modes in the CC case, for several values of 
the asymmetry parameter t. The low density solution represents the unstable isoscalar 
mode (black line). The remaining stable modes are depicted with red or blue lines for 
isovector or isoscalar character, respectively. 



